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Abstract 

We analyze K — K° oscillations in space in terms of propagating wave packets 
with coherent K$ and Kl components. The oscillation probabilities Pk°-+K°{ x ) an d 
Pk°^k°( x ) depending only on the distance x, are defined through the time integration 
of a current density j(x,t) . The definition is such that it coincides with the experi- 
mental setting, thus avoiding some ambiguities and clarifying some controversies that 
have been discussed recently. 



1 Introduction 



Due to ongoing and planned experiments in kaon physics there has been a renewed 
interest in space-time oscillations of neutral meson systems like K°-K°. In particular, 
in ||, |H an effort has been made to compute oscillation probabilities, like PK°-+K°( r ) 
and P K o^xo(r), which depend now solely on the spatial distance r between the K° 
production point and the position where K° or K° are detected. The expression P(r) 
mimics then closely the usual experimental circumstances where distances rather than 
times are measured. This issue can be studied in single neutral kaon production in 
reactions like ir~p — > AK°, K + n — > pK°, pp — > K + ir~K°, or in kaon pair production 
in reactions such as e + e~ — ► KgK^ at factories like Daphne ||, ||. We will 
concentrate in the present paper on single kaon production — a problem which finds 
its analogy in recent discussions on neutrino oscillations @ 3 §] ■ - and point out the 
differences to a ^-factory at an appropriate place. 

The methods and arguments given in [jl], [2], |j for the single kaon oscillations (for 
the analogous neutrino case, see @, to justify a particular derivation of the spatial 
probability P(r) are rather different from each other and so are partially the results. In 
our opinion, the origin of the discrepancies (in both, the derivation and the results) can 
be traced back to the absence of a clean quantum mechanical definition of a probability 
which depends only on the distance r. If the question of what is the probability to detect 
either a K° or a K° at a distance r from the kaon production point is a legitimate one, 
then a quantum mechanical definition should indeed exist. Sometimes the classical 
formula r = vt is invoked to transform P(t) into P(r). However, strictly speaking, this 
procedure might not be appropriate in a quantum mechanical context due to the finite 
spread of the wave packets. Furthermore, in the neutral kaon system the existence of 
a different mass for each component Kl and K$ implies different mean velocities. 

From the point of view of quantum mechanics a probability has to be defined in 
terms of the wave function ty{x,t) which encodes all the information of the system. This 
probability should then be constructed according to our experimental requirements. A 
position dependent probability P(r) for propagating states can only be obtained by 
integrating an appropriately defined quantity over the time variable. Following the 
above arguments one concludes that the quantity to be constructed in the first place in 
terms of *&(x,t) should be a probability density in time (dependent on x and t) which 
corresponds to the current density. Indeed a physical detector at a distance r from 
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the production point measures the time integrated flux of probability flowing across its 
surface , i.e., 

P(r)= f dA f dtj(x,t) (1) 



A 
with 

where A is the surface and n is a unit vector pointing outwards. 

The problem then is to find a definition for the current in the case of mixed states 
such as K° and K°. Note that in quantum mechanics the current j(x,t) is usually de- 
fined only for states of definite mass, nevertheless we will show that the corresponding 
generalization is possible in our case. In doing so we will restrict ourselves for sim- 
plicity to 1 + 1 dimensions. Moreover, since the controversial point in connection with 
PK°-*K°( r ) aim Pga^Koir) refers to neutral kaon propagation in space-time and it has 
neither to do with relativistic effects nor with CP-violation in the neutral kaon system, 
we can work in i) a non relativistic limit for which the current take simple expressions, 
it) a CP conserving theory, and Hi) a situation with stable kaons, Tg = Tl = 0. The 
last two points are consistent with the Bell-Steinberger relation in the limiting case of 
(Ks\Kl) = |J. Note that the non-relativistic limit refers to retaining the rest mass 
in the energy-momentum relation, while neglecting for the moment higher orders of 
v/c. The modified Schrodinger equation is then: 

i§;*s/L(x,t) = ^—^^ s/L (x,t)+m s/L ^ s/L (x,t), (3) 

which is the appropriate non-relativistic limit of either Klein-Gordon or Dirac equations. 
The relativistic counterparts including finite width effects iXs/L 7^ 0) wm be easily 
inferred from our results. 

In next section we recall some basic formulae of the neutral kaon states and of 
wave packets. In section 3 the current density is obtained and in section 4 we compute 
the oscillation probabilities. Some brief conclusions follow. 



2 Wave packets for K° and K° 

In the neutral kaon system one can distinguish two sets of states. The "strong- 
interaction" basis is given by the states | K°) and | K°), with well defined strangeness 
but undefined masses, fulfilling the orthogonality condition (K° \ K°) = 0. In contrast, 
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the "free-space" propagating states | Kg) and | Kl) have well defined masses, m s / L , 
but are not strangeness eigenstates. This is the appropriate set to describe time evolu- 
tion in free space given by \K S/L (t)) = e~ lXs J Lt \K s/L ) with X s/L = m s/L - i/2F s/L -> 
ms/L m our limit. The relation between both sets is given by 

\K s/L )=p\K°)±q\K0) (4) 

with | p | 2 + | q | 2 = 1 and | p | 2 — | q | 2 ^ if CP is violated in the mixing. 
A general space-time dependent state \K S / L (x,t)) can be written as 

\K s/L (x,t)) = 9 s/L (x,t)\K s/L ). (5) 

From Eqs. one can now obtain the time evolution of an initially produced K°. 

In the two dimensional space spanned by | K°) = ( ] and | K°) = I I one has 



K°(x,t)) 



<S> K o(x,t) \ _ 1/ <St s (x,t)+y L (x,t) 
**oOM) J 2 V l(*s{x,t)-*L{x,t)) 



(6) 



$ S (x,t)-y L (x,t) 

where the last term in (|6|) holds only in the CP conserving limit p = q = 1/ \/2. In case 

(^5 — ^>L 
-^S-^L 

The wave functions ^^/^(x, t) are taken to be wave packets which are known to 
give a realistic description of the space-time evolution. We choose for simplicity a 
Gaussian form 

*S/L(*,t) = 7^77 Pjke-^-^^e^s/dm^ (7) 
where according to (||) 

k 2 

E S /L = u S / L (k) = m s/L + , (8) 

Zm S/L 

and kg/L are the mean values of the momenta of the short/long kaon states. Recall 
also that a is related to the width of the packet, Ax = a/2, and that the uncertainty 
relation gives AkAx = 1/2, i.e., Ak = 1/a for both K$ and Kl- Upon integration over 
k an explicit formula for ^g/L{x,t) can be obtained 

/ 2 \V4 1 f (x-v s/L tf 1 . 
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with 

Vi = -Oi + hx - LOi{ki)t + W7- - i = S,L (10) 

u(k) defined in Eq. (@), and 

2 k- 

li = ~, Vi = — , tan(26»i) = 7^. (11) 

mjcr mi 

Eqs. ®-([[l]) represent the usual Gaussian wave packet displaying a well-known spread 
in time for definite mass eigenstates K$ or K^. 

Note that in our single K° production processes, energy-momentum conservation 
requires that the two Kg/^ components of the K° wave packet should have slightly 
different mean energies Eg /l and momenta k^jL due to the mass difference 5m = 
mi ~ m s > 0- in this sense we fully agree with the main point raised by Srivastava 
et al.[j^], also advocated in |7j ; concerning the presence of two different energies and 
momenta. Both mass eigenstate components have to be produced in the coherent 
superposition l/2{^>s( x , t = 0) + ^(a^i = 0)] corresponding to an initial K°. This 
coherent production of overlapping wave packets of K$ and also requires 

Ak»5k = k L -k s (12) 

In this respect the situation is very similar to that considered in neutrino oscillations 
whose wave packet treatment has been discussed in || . It is also interesting to observe 
that in our treatment single K° production somehow mimics a two-slit experiment in 
momentum space, 5m playing the role of the "separation" between the Kg and Kl 
"slits". 

There is a subtle point worth mentioning here referring to the Bargmann super- 



selection rule [[R|,|ll]. This rule states that in non-relativistic quantum mechanics 
Galilean invariance forbids a coherent superposition of states with different masses. 
We do not consider this a serious problem here. Beyond the non-relativistic limit such 
restrictions do not hold as is clearly experienced in the neutral kaon system. Note that 
Eq. (^j) is a limit of a relativistic equation and that we retain the different rest masses. 
A relativistic current has then to reduce to its non-relativistic counterpart which is 
based on Eq. (||). Notice also that Galilean invariance cannot be strictly maintained 
if, for instance, we would use time dependent potentials in the Schrodinger equation 
(like space-time dependent electromagnetic potentials) [|ll]]. 
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3 Currents for K° and K° 



The current density for the K^/l niass eigenstates whose wave functions obey Eq. (|3|) 
is defined as 

is/L&t) = (% /L (x,t)^ s/L (x,t)) , (13) 



??? 



S/L 



and satisfies the usual continuity equation j^ps/L + -§^3s/L = 0> with Ps/l{ x ^) = 
^S/l^s/l- For the K° and K° states we expect that the corresponding currents are 
not conserved as these states mix in time due to the mass difference 5m in their K$ 
and Kl components. Thus, the modified continuity equation should read 

d d , . . 

qIPk°,K° + q^3ro,KO - d K°,K°> V 14 ) 

where p K o = ty* KQ ty K o, p^o = ^*^o^k° and d R0 is a kind of "diffusion term" to be 
determined. 

Let us list some general requirements for the current density j K0 R0 '- i) It should 
contain only "velocity terms", i.e., terms of the form l/rrii y&j^'&k) , with k = S,L, 
whereas one has no similar condition for d K0 R0 . ii) As the mass difference goes to zero, 
5m = m^ — ms — ► 0, we should have d K0 R0 = 0, independently of the form of the wave 
packets ^ s/L- l n this limit obviously one has j K o = = js = 3L- I n the zero 
decay width limit the total number of neutral kaons must be a conserved quantity, i.e., 

j f (Pro + P K o) + ^ {j K o + 3 K o) = (15) 

Hence, d K o = —d^ should hold. 

We can now derive the expression for the current density j K0 k°- First, we take 
the time derivative of the densities \^k° k°^k° k )' where the wave functions are 
given m Eq. d). We then use the Schrodinger equation ^) and separate the "velocity 
terms" according to point i) above. The linear combination of such terms will be called 
j K o ko ■ All other terms that are not of velocity type are included in d R0 R0 . We thus 
obtain 

Jko,ko = 1 Us + 3L) + ^JkIko, (16) 

where j s / L are the usual current densities (|13|) corresponding to the Kg/ L states, and 
?t?n r ?n contains the interference terms 



■int _ 1 

J K o ~ 2 



±lm [n^ L ) + ±lm (n&s)] = d7) 
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For the "diffusion terms" we get 

d K0 = J (-L - -L) M^ S ) + - roj? )/m(n*x) = -^o (18) 

It is now easy to check that the three aforementioned requirements on the current 
densities are indeed fulfilled. In particular, d K0 R0 — > as 5m — > 0, and Eq. (15) is also 
trivially satisfied due to the antisymmetry properties of j % J% ^ and d R0 R0 in (17) and 
(18), respectively. 



The current densities fllq ) can also be written in terms of the wave functions 

and 

j& = \ ~ 2 _\ 6m2 (2mlm^±^) - SmIm^ ^ I<0 )) , (19) 

where m is the mean value of the mass, m = (ms+mt)/2. Note that in contrast to the 
densities p R0 each of the two currents j R0 R o contains both wave functions x°- 
This means that the constructed currents are unique in the sense that all "velocity 
terms" are contained in j K0 R o and not in d R0 R0 . 

4 Distance Dependent Probabilities P^^k K°( r ) 

Given the currents j R0 R0 in ( |l6| ) the x-dependent probabilities 

P K o^ K o(x) = / dtj K o(x,t), 
Jo 

POO 

P k o^ko(x) = / dbj&fat) (20) 
JO 

can now be calculated in analogy to Eq. ([!]) (recall that for simplicity we are working in 



1+1 dimensions). Before evaluating the integrals of Eq. ( P0|) it is convenient to elaborate 
on the kinematics involved. Note that we will use consistently the non-relativistic limit, 
Eq. (||), and restrict ourselves to the case of single kaon production in "two-to-two" 
reactions like tt~~p — ► AK°, K + n — > K°p, etc. Working in the CM-system with the 
total energy y/s = m + M + Q, where M is the final baryon mass, one has 

5k = k L -k s , k = -{k s + k L ) 

kSk ( k 2 \ 1 

5 E = E L -E s ~— + 6m[l-^}, E^-(E L + E S ) 
m \ Im* I 2 




1 

m z 2 



Sv = v L -vs~ — - 2 -, v = -(v s + v L ) (21) 
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The Q-value of the reaction can be chosen in the very wide range 



m » Q » 5m 



(22) 



imposed by our non-relativistic treatment and the need to produce a K° in a coherent 
Kg/L superposition. Expanding the relevant variables one finds 

2 



A: 



E 



2Mm 
M + m 



Q, k5k 



-dm 



Mm 
M + m 



M 



M + m 



m + 



m 



M + m 

and, more importantly, 

5v 5m 
v m 



Q, 5E ~ 5m 



M 



Q 



M 



M + m (M + m) 2 



Q 



M + 2m m 
M + m + Q 



O 



5m\ 

~q) 



o 



5m 



mv 



(23) 



(24) 



We stress that in spite of working in the non-relativistic limit v cannot be arbitrarily 
small as to imply a Q-value (or a kaon kinetic energy) of the order of the mass difference 
5m. In such unrealistic situation one would be favouring the production of the Kg 
component over the slighlty heavier Kl one. Taking into account that 5m/m ~ 10~ 15 
Eqs. (22) and (p4j) imply v » 1CP 7 . Notice, however, that even velocities of the order 
of 1CP 3 would require an unrealistic high degree of fine tuning of the CMS energy y/s 
(Q ~ keV). 

At this stage it is instructive to point out some differences to a ^-factory like 
Daphne. In the latter case one usually studies correlation probabilities of K$ and Kl 
decays depending on two space-time points ||. The kinematics dictates in this case 
5v/v = —(5m/m) rather than (|24j). 

In terms of the Gaussian wave packets 
be cast into a simple form 



Jt 



Vi + xjft 



7T o(l + Tft 2 ) 3 / 2 

Introducing the new variables 



exp 



, the currents jg/L appearing in ( J16| ) can 

(X - Vit) 2 



'a 2 {l + 1 2 t 2 ) 



Vit 



L,S 



(1 + 7^ 2 ) 1/2 
we can compute to a very good approximation 



dtji(x,t) 



»i 



iafci<<-l 



dzi e 2z ? ~ 1 



(25) 



(26) 



(27) 



8 



Indeed in the upper limit x » a since the position x of the detector is considered 
to be far from the K° production point, far in comparison to the width of the packet 
a/2. Otherwise, one would be measuring the inner wave packet effects rather than its 
space-time propagation which is our interest here. In the lower limit of the integral 
the momentum k is much larger than its dispersion Ak = 1/a for any realistic velocity 
v ~ 1CT 3 - HT 1 . 

Neglecting in ( |l~7|) polynomial terms of order 5m /m and expanding further the 
remaining expression in leading order of 5m/ m, we can write 



■int 



v + X7 t 



-exp ■ 



vtf 



t 2 5v 2 



vra(l+ 7 2 t 2 ) 3 / 2 
cos 9{x,t) + sin 6(x, t) 



a 2 (l + 7 2 t 2 ) 
Swyt 



exp ■ 



2a 2 (l + 7 2 t 2 ) 



where 



6(x,t) 



-5Et + 5kx + 



2(v + x-/ 2 t)_ 

2t5v(vt - x) 



(1 + 7 2 * 2 ) 



(28) 



(29) 



is the crucial phase and the mean 7 = 2 /ma 2 , v = k/m have been defined in analogy 
to Eq. (|TT|) . Observing the similarities between the Eqs. (p5[) and (|2§|), we introduce a 
variable as in (pq), az = ^^2^2*1/2 1 whose inverse in the region of interest < t < 00 

- azy^v' 2 + 7 2 (x 2 — a' 2 z 2 ) 



is 



t(z) 



(30) 



(v 2 — 7 2 a 2 ^ 2 ) 

for f > and x > 0. The integrand has now a form of e~ 2%2 f(z) which can be expanded 
in a Taylor series around z = 0. Note that i(0) = x/v and that the coefficient of sin# 
in (28) is bounded by 1/2 (5v /v) 1+ ffiffi/ ) — 1/4(5^/^) ~ 0(5m/Q) and hence small by 
virtue of Eq. (|22|). In leading order we then obtain 



dtf$(x,t) ~ 



cos 0(x, t(0))exj» 



t 2 (0)ow 2 



2a 2 (l + 7 2 t 2 (0))w 2 



~ cos9(x,x/v) 



(31) 



2 2/2 

since the additional Gaussian factor has a negligible exponent, l/2(5v' 2 /~/' 2 ) i+^zjy v s < 
l/2(5u 2 /7 2 ) ~ l/2(5k/Ak) « 1 due to @. 

Taking into account the kinematics given in Eqs. (^2| 

be 



8(x, x/v) 



-5m I 1 



1 



the phase turn out to 
(32) 
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Putting now everything together we obtain 

P K o^ K o t Ko(x) = - (l ± cos 



5m— (1 v 2 ' 

v 2 ' 



(33) 



which is our final result in non-relativistic approximation and for stable kaons. Rein- 
troducing the non-vanishing T s /l decay widths and higher orders in v/c, Eq. (|33|) can 
be cast into its final form as 



+ e 



db e 2 



^(rs+r L ) 



cos [5mr 



(34) 



which corresponds to the usual expression in terms of the proper time r. 



5 Conclusions 

We have argued that in order to avoid theoretical ambiguities it is essential to use a 
proper definition of space-dependent probability P(r): the time and surface integrated 
current density, j(x,t). We have shown that it is possible to construct such a current 
density for coherent superposition of different mass eigenstates. Further, we have used 
throughout wave packets which describe the physical space-time evolution and take 
into account the necessary overlapping of wave functions to obtain interference effects. 
Wave packets have also been used in Q to compute space-dependent decay correlation 
probabilities at a ^-factory. Their results are quite in line with those obtained in the 
present paper. 

In our treatment the two components K$ and Ki of the neutral kaon evolve 
in space-time with different momenta and energies as dictated by energy-momentum 
conservation in a situation that can be visualized as a two slit experiment in momentum 
space. We have then obtained the oscillatory term cos [<5mr] , which corresponds to the 
usual expression but without the recourse to any classical formula to transform time 
probabilities into space probabilities. 
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